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Abstract 

We reconsider Ising spins in a Gaussian random field within the 
replica formalism. The corresponding continuum model involves sev- 
eral coupling constants beyond the single one which was considered in 
the standard 4 theory approach. These terms involve more than one 
replica, and therefore in a mean field theory they do not contribute to 
the zero-replica limit. However the fluctuations involving those extra 
terms are singular on the Curie line below eight dimensions, and by 
the time one reaches the dimension six, it is necessary to keep them in 
the renormalization group analysis. As a result it is found that there 
is no stable fixed point of order (6 — d). Whether this means that 
there is no expansion in powers of (6 — d), or that the transition is 
driven to first order by these fluctuations, is difficult to decide at this 
level, but it explains the failure of the (d, d — 2) correspondence. 
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In spite of more than twenty years of active and continuous research, the 
simplest disordered system, the random field Ising model (RFIM), is still far 
from a definite understanding. The initial (d, d-2) correspondence, between 
the RFIM in dimension d and the pure system in dimension (d-2)(jl], |2], [3|]), 
was suspected from the very beginning to break down in lower dimensions, 
if one followed the arguments of Imry and Ma ([|J). The correspondence 
was definitely killed when the mathematical approach of Imbrie, and sub- 
sequently Bricmont and Kupiainen, finally established the existence of an 
ordered state at low temperature in three dimensions ([|5], §). But this still 
left room for various possibilities. The elegant "derivation" of the (d, d- 
2) relation by Parisi and Sourlas (||) made it clear that non-perturbative 
issues could destroy the essential underlying supersymmmetry (since a Jaco- 
bian, the absolute value of a determinant, was replaced by the determinant 
itself which could change sign). One could thus imagine that the (d, d-2) 
correspondence for the critical exponents, although true to all orders in an 
expansion in powers of e = 6 — d, was modified by terms, such as exp — 1/e, 
which might very well turn out to be quite significant at e = 3. However 
progressive evidence was provided that the situation could be more serious 
than that. Numerical experiments in the supposedly paramagnetic phase, 
pointing at large time scales before one reached the ordering temperature, 
as well as the theoretical approach of Mezard and Young (0), led to the 
idea that the phase diagram itself had to be reconsidered. Coming from the 
high temperature regime it seems that one meets a glassy phase, marked by 
replica symmetry breaking (RSB), in a region in which there is no magnetic 
ordering, as in a spin glass. We shall present in a subsequent publication 
arguments which substantiate this analysis. However here we argue that 
within the original models of Ising spins, or of the 4 theory, in a Gaussian 
random field, another scenario might take place. Indeed we show, within the 
replica approach, that new interaction terms should be considered. These 
terms naively disappear in the zero replica limit, but it is argued that below 
dimension eight, these new interactions develop singularities near the Curie 
line in the n = limit, which change the analysis when one reaches the di- 
mension six. The critical renormalization group flow, if performed as usual 
on the Curie line, contains then five coupling constants, and it is easy to show 
that there is no stable fixed point of order (6-d). Above the Curie line these 
singularities are cut-off by the finiteness of the correlation length, and there 
is a non-uniformity in the approach of the Curie line, versus the zero-replica 
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limit. We are not able to say whether the absence of a stable fixed-point 
means that there is a non-perturbative fixed point elsewhere or whether the 
transition becomes first order, but it shows that the dimensional reduction 
may break down even at lowest order. 

1 The basic action 

We consider successively, within the replica formalism, 4 and Ising spins in 
an external Gaussian random field h(x) with 

< h(x) >= 0; < h{x)h{y) >= A5(x - y). (1.1) 

Therefore starting with 

S(<f>) = J d^(V0) 2 + ir o 2 + |0 4 - h(x)<f>(x)} (1.2) 

(the Boltzmann weight is exp — S) and replicating n times we average over 
the random field and obtain 

r n 1 1 n A 

S{<t> a ) = J <M£(-(V0 Q ) 2 + -r^l + Ul) - -Y,M*mx)] (1-3) 

a=l " a,fS 

Let us introduce the notation 

n 

a k = J2(Ux)) k ; (1.4) 

a=l 

the symmetry of this action under permutation of the replicas allows for five, 
a priori relevant, coupling constants corresponding to the quartic interactions 

04, CT10-3, (or 2 ) 2 , (0l) 2 02, (0l) 4 • 

Indeed if we were considering the renormalization programme for the 
various correlation functions of the theory at finite n , we would be forced to 
introduce these coupling constants, even if they are not all present initially. 
In order to substantiate this claim, let us start with Ising spins and derive 
the corresponding c/> 4 theory in the usual way. We begin with 

(3H=-W JiASj - E kSi, (1.5) 
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and use the Hubbard-Stratonovich identity 

E JiASj+Y, h ^ = c j (IW ex p[-^ E &^V;+E(&+ Wl 

(1.6) 

We then introduce n replicas and sum over the spins. This gives the partition 
function 

. N n i 

Z" = c/nn^exp[--^^ a J^,a + E lncos M4,a + ^)] (1-7) 

i=l a=l ij,a i,a 

We then expand in powers of , up to degree four: 

ln £2 ^ 1 - ^ + t * + V + ¥ + ^ + °^ ^ 

in which r = tanh h, and average over the random field h: 

cosh(0 a + h) .1. . r 2 2 

< exp E In ^ >= exp(-(l - r> 2 + -a\ 

- — (1 - 4r 2 + 3r 4 y 4 + ^(r 4 - 3r 2 2 )^ + -(r 4 - r 2 > 2 2 

+ i(r 4 - r 2 )<Ti<T 3 - i(r 4 - r 2 >?(7 2 + O(0 6 )) (1.9) 

with 

r p = / d/wT^ (cosh /i) n (tanh /i) p (1.10) 

v27rA J-oo 

At the level of the quadratic terms (1— r 2 )<7 2 shifts the constant interaction 
term J~ l {q = 0) and gives an effective 'mass'. The other term t 2 J]0o0/3 is 
due to the average over the random field. For a small random field, i.e. small 
A , in the n = limit 

r 2 = A - 2A 2 + 0(A 3 ); r 4 = 3A 2 - 20A 3 + 0(A 4 ) (1.11) 

We are then led to a 4 - theory with five coupling constants: 

(3H = /^(lE[(v0.) 2 + ^]-f E< 



■^j- CT 4 + ^-aia 3 + ya 2 + — ^^2 + — aj (1.12) 
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A priori the four additional terms beyond the single interaction <j\ , which 
is the only one usually present, involve more than one single replica. There- 
fore if one writes the quenched average of the replicated partition function 
in terms of W = LogZ , 

< Z n >= exp[n <W> + — (< W 2 > - < W > 2 ) + 0(n 3 )], (1.13) 

it is easily seen that the contribution of the four additional coupling constants 
are of relative order n (or higher) and thus could be discarded. However this 
argument is true provided the coupling constants Ui for i > 1, are not singular 
in the n = limit. 



2 Singularities below dimension eight 

At the level of mean field theory the coupling constants that we have obtained 
in (|1.9| ) are all finite in the n = limit. However the situation changes 
with the fluctuations: consider the one-loop fluctuations in the theory at the 
critical point. These fluctuations are given by the determinant generated by 
the Gaussian fluctuations around the mean field, namely 

exp-iTrln(l + G S {2) ) (2.1) 



in which 



Gf(q) - {q ^ t) + {q 2 + t)iq ^ +t _ nA) ( 2 - 2 ) 



and is the second derivative of the action ( |1 . 1 2 ) with respect to (fi a and 



Sap = M^l^a + U 2 C?10 a + -U 3 a 2 + ^U 4 (Ti 2 ] + ^U 2 (4>1 + ^) 

+U 3 <j) a <f) l3 + ^M 4 CT 2 + U 4 CTi(0 a + <pp) + ^U 5 0"i 2 (2.3) 



in which we have used the notation (|1.4[) . The fluctuations terms which are 
quartic in <p are obtained by expanding ( |2.1| ) to second order in S^ 2 \ We first 
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consider the most singular part in which for the propagator Go we retain 
twice the A part of ( |2.2|) . It involves the sum 

A 2 (E S S) 2 = A 2 [^i + 4/m 2 + nu 3 + n 2 u 4 )a 2 

5 

+ O2 + W3 + -nu 4 + n 2 u 5 )a 2 } 2 (2.4) 
multiplied by the one-loop integral 

1 



i n (p,t) = J A- 



4(0, nA) = / d d q ~ (nA)-^r. (2.6) 



(g 2 + t)(g 2 + t - nA)((p - g) 2 + i)((p - g) 2 + t - nA) ' 

(2.5) 

The contribution of these fluctuations to the effective coupling constants 
are thus related to I n (0, t). Note that one cannot let t go to zero first, since 
this would yield a pole of the integrand on the line q 2 = nA. Therefore in 
order to determine the critical behaviour of this effective interaction we have 
to examine I n (0, t) for t of order nA, and then let n go to zero. On the line 
t = nA, which is close to the usual critical t — line, the integral l n (0,t) 
becomes singular, in dimensions lower than eight. A simple power counting 
gives 

1 / S-d 

q 4 (q 2 + nAy 

Therefore these fluctuations yield singular contribution to the coupling con- 
stants contained in ( |2.4|) in dimensions lower than eight in the zero-replica 
limit. Thus the coupling constant M3 receives a contribution proportional to 
~u\ in dimension six, likewise for u 2 - On the other hand M4 behaves like 
Uiu 3 /n, i.e. like 1/n 2 and u 5 like u\jn, i.e. 1/n 3 . Consequently it is not 
legitimate to discard the coupling constants which couple several replicas in 
dimensions smaller than eight. 

It is interesting to analyze how a dynamical approach, based for instance 
on a Langevin formalism, would recover those singularities, despite the ab- 
sence of replica. In the Martin- Rose- Siggia formulation of the dynamics (H), 
there are two fields, namely 0(x, t) the order parameter, and the conjugate 
field <f)(x,t). The equilibrium 4 theory leads to a single dynamical interac- 
tion term of the form 00 3 . However fluctuations build up additional couplings 
such as 2 2 ,etc, with singular coefficients in dimensions lower than eight. 
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It is then necessary to introduce a waiting time and invariance under time- 
translation invariance is broken . The n = singularity in the critical domain 
is then replaced by a singularity when the waiting time becomes large. A 
full comparison between the dynamical approach and the equilibrium replica 
formalism goes beyond the scope of this letter. 



3 The renormalization group flow near six di- 
mensions 

Since the multi-replica coupling constants are singular below dimension eight, 
one cannot discard the new interaction terms in the zero- replica limit, since 
they blow up as powers of 1/n in dimension six. Therefore we keep now 
the full action ( 1.12D with the five coupling constants, and repeat the usual 



renormalization group approach near six dimensions. We know that if we 
had kept only the coupling constant u\ we would have found a stable fixed 
point of order e = 6 — d, which would yield a critical behaviour identical to 
that of the pure system near dimension four(|| [J). At one-loop we obtain 
the renormalization group equations by returning to ( |2.1| ), expanding again 
to second order in S^ 2 \ but now keeping only one A in the propagators. The 
result is proportional to 

XXSI S apT = + nu 2 ) 2 a 4 + (ui + nu 2 )(u 2 + u 3 + nu^a^ + 



a p 



+ -^[n(u 2 + u 3 + nu 4 ) 2 + 2(u 2 + u 3 + nu±){u\ + nu 2 )\o 2 

+ [2(u 2 + u 3 + nu A ) 2 + ^(3m 4 + nu 5 )(ui + 2nu 2 + nu 3 + r?u^)\oio 2 



+ [{u 2 + u 3 + nM 4 )(3M 4 + nu 5 ) + -(w 4 + nu 5 ) }o l 



multiplied by an integral proportional to 1/e. Redefining the coupling con- 
stants to contain geometric factor ( surface of the unit sphere divided by 
(2n) d ) we obtain the five one-loop beta functions : 



(3.1] 



(3i = —eui + 3A(w! + nu 2 ) 2 

(3 2 = -eu 2 + 3A(«! + nu 2 )(u 2 + u 3 + nu 4 ) 
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(3 3 = -eu 3 + A[n(u 2 + u 3 + nu 4 ) 2 + 2(u 2 + u 3 + nu 4 )(u 1 + nu 2 )} 
(3 4 = -eu 4 + A[4(w 2 + u 3 + nu 4 ) 2 

+n(u 2 + u 3 + nu 4 )(3u 4 + nu 5 ) + (3u 4 + nu$){ui + nu 2 )} 

(3 5 = - eu 5 + 12 A[(u 2 + u 3 + nu 4 ) (3u 4 + nu 5 ) + -(3u 4 + nu 5 ) 2 ] (3.2) 

These equations make it clear that the renormalization mixes the five 
operators that we are considering. If, notwithstanding the potential singu- 
larities of u 2 , ■ ■ ■ ,w 5 , one lets n go to zero first in the equations |T2|, one 
recovers the usual dimensional reduction fixed point ^1 = 3^ + 0(e 2 ),u 2 = 
■ • ■ — «5 = 0. It is immediate to check that this fixed point is unstable and 
that the u^s for i > 2 grow. Furthermore we have seen in the previous sec- 
tion that the effective interactions with the symmetry of those extra coupling 
constants, are indeed singular when n goes to zero. Therefore instead of let- 
ting n go to zero first, we give to these coupling constants the dependence 
in n inherited from these singularities, i.e. we combine the n-dependence 
from these singularities with the sum over replica indices. Examining those 
beta functions ft.2|, it is immediate to verify that if we redefine the coupling 
constants 

#1 = Ami, 92 = nAu 2 , g 3 = nAu 3 , g 4 = n 2 Au 4 , g 5 = n 3 Au 5 (3.3) 

the number of replicas n drops completely from the equations. Note that 
the dependence in n of these coupling constants is exactly consistent with 
the analysis of the effect in dimension six of the multi-replica singularities 
analyzed in the previous section. This leads to the important conclusion 
that the existence and stability of the fixed points is now independent of the 
number of replicas. In terms of the couplings gi we now have 

Pi = -egi + 3(5-1 + g 2 ) 2 

= -£Q2 + 3(^1 + g 2 ) (g 2 + g 3 + g 4 ) 

fa = -e# 3 + (92 + 93 + 9a) 2 + 2(^2 + #3 + 04)(0i + 92) 

(3 A = -eg 4 + A(g 2 + g 3 + g 4 ) 2 + (g 2 + g 3 + g 4 ) {3g 4 + g s ) + (3g 4 + g 5 ) (g 4 + g 2 ) 

P5 = -e#5 + 12((?2 + 9s + 94)(Sg 4 + g 5 ) + 3(3g 4 + g 5 ) 2 (3.4) 

(we have kept the same name to the beta functions, in spite of an obvious 
refinition by a multiplicative factor corresponding to the redefinition fl3.3| )). 
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We now look for fixed points gi* which are zeroes of those equations and 
satisfy the stability condition given by the non-negativity of the eigenvalues 
(or rather the real part of the eigenvalues if they were complex) of the matrix 

fijj — 7: -|g* (3-5) 

dgj 

at the fixed point. Note that if g\ was the only coupling, the fixed point 
gi* = | would be stable, but it is no longer true when the others are present 
since ^ = — | being negative, the matrix cannot be non- negative. The 
discussion of the fixed point solutions is easier in terms of the variables 

g\+92 = x, g 2 + g 3 + #4 = y, 3# 4 + g 5 = z. (3.6) 

and the first equation ( |3.4| ) leaves us with two possibilities a) x = and 
then ^ is negative, b) x ^ , then x + y — | and |^ = — e + 2(x + y) 
is negative. We then conclude again that there is no stable fixed point. 
Therefore, whether we take into account the singular dependence in n of 
the effective interactions or not, one does not find a stable fixed point of 
order (6 — d). This still leaves us with two possibilities (i) a stable fixed 
point which is not of order e (ii) a runaway flow indicative of a first order 
transition. However in both cases there is no epsilon-expansion of the critical 
properties of the random field Ising model, and the dimensional reduction 
breaks down at first order in e. 



4 Branched polymers 

One may wonder at this stage whether the previous analysis would not en- 
danger as well the beautiful results of Parisi-Sourlas on dimensional reduction 
for branched polymers(||), for which there are only reasons to believe that 
it works. We shall not repeat here their original derivation but simply take 
the replicated field theory version of the problem, which is now a 3 theory. 
The action is similar to ( |1 . 1 2[) 

(3H = /^B(V0.) 2 +^]-f E ^+|^3+f^ 2 +|k 3 ). (4.1) 

Again the theory has been based on the single coupling constant u\ since the 
other two drop from the problem in the zero- replica limit. The upper critical 
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dimension is then eight, and the authors of (0) proved that the epsilon- 
expansion about eight dimensions (based on a purely imaginary coupling 
constant) is identical to the (6 — d) expansion for the 3 theory. 

Similarly here the multi-replica coupling constants develop singularities 
in the zero-replica limit below dimension twelve which make it necessary to 
reconsider the analysis of the (8 — <i)-expansion. Following the same lines we 
have obtained the three beta functions. Again n drops from the determina- 
tion fo the fixed points if we redefine 

#i = Au u g 2 = nAu 2 , g 3 = n 2 Au 3 (4.2) 

and we find 

Pi = ~9x - 2 9i(gi + 92) 2 - ^gi(gi + 3# 2 + gs) 2 

P2 = ~^92 ~ ^92(9i + 3# 2 + 0s) 2 

-^9i(9i + 92)(2g 2 + g 3 ) - 2g 2 (g x + 02X01 + 3g 2 + g 3 ) 
As = —zga + 93(91 + 92) 2 - 3(25-2 + 03) 2 (fi'i + 92) - 602(01 + 3# 2 + g 3 )(2g 2 + g 3 ) 
13 

-^93(91 + 3fi- 2 + 03) 2 - 602(01 + 02) (202 + 03) (4.3) 

The Par isi- Sour las fixed point is g±* 2 = —if, 02* = 03* = 0, and it is im- 
mediate to verify that in this case it is (marginally) stable. Therefore the 
mechanism which invalidates the epsilon-expansion for the RFIM is harmless 
for the problem of branched polymers. 
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